Lyapunov exponent for the laser speckle potential: a weak disorder expansion 
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Anderson localization of matter waves was recently observed with cold atoms in a weak ID 
disorder realized with laser speckle potential @. The latter is special in that it does not have 
spatial frequency components above certain cutoff q c . As a result, the Lyapunov exponent (LE), 
or inverse localization length, vanishes in Born approximation for particle wavevector k > iq c , 
and higher orders become essential. These terms, up to the order four, are calculated analytically 
and compared with numerical simulations. For very weak disorder, LE exhibits a sharp drop at k 
= ig c . For moderate disorder (a) the drop is less dramatic than expected from the fourth order 
approximation and (b) LE becomes very sensitive to the sign of the disorder skewness (which can be 
controlled in cold atom experiments). Both observations are related to the strongly non-Gaussian 
character of the speckle intensity. 

PACS numbers: 42.25Dd, 03.75.Nt, 72.15.Rn, 05.10.Gg 



Technological progress in experiments with ultracold 
atoms provides an extraordinary level of control, allow- 
ing investigation of various quantum phenomena [l| . One 
such phenomenon is Anderson localization (AL) || of 
matter waves. Formulated originally to explain the ab- 
sence of spreading of quantum-mechanical wave function 
in a disordered potential, it was recognized later as a 
common signature of wave propagation in random media, 
where waves may become exponentially localized because 
of the destructive interference between its multiple scat- 
tered components. Some indications of AL were observed 
with light [|[ , microwaves [H and ultrasound @ , while in 
electronic systems it is hindered by finite temperature 
dephasing and interactions. 

Cold atoms offer a unique possibility to study matter 
wave localization in conditions of very low temperature 
and tunable interactions (not available in solids) , absence 
of absorption and controllable disorder. The latter may 
be introduced by different techniques, one of which is us- 
ing static laser speckle, whereas potential felt by atoms is 
proportional to the speckle intensity with the sign of the 
detuning from the atomic transition [y]. Laser speckle, 
generated by passing expanded laser beam through diffu- 
sive plates, are special in that they have (i) exponential, 
i.e. strongly non-Gaussian, intensity distribution and (ii) 
finite support of their power spectrum Q ■ As discussed 
below, these properties have strong effect on the local- 
ization properties in fD systems. 

Recently, AL was observed in Bose-Einstein conden- 
sates expanding in weak ID disorder realized with laser 
speckle 8] . A simplified model for such experiments con- 
siders expansion of the condensate released from a har- 
monic trap as occurring in two stages (1) explo- 
sive conversion of the interaction energy into kinetic en- 
ergy, in which weak disorder may be neglected, followed 
by (2) expansion of the non-interacting gas in the disor- 
dered potential. The first stage ends up with a certain 
momentum distribution of the atomic cloud (TlT | , which 
provides an initial condition for the second stage. Thus, 
the problem reduces to basically a single-particle local- 
ization, characterized in ID by Lyapunov exponent (LE), 



or inverse localization length [12|. The latter, in a weak 
f D speckle disorder with a spatial frequency cutoff q c , has 
some peculiarities, because in ID (as opposed to higher 
dimensions), elastic scattering is a "binary" process: ei- 
ther the particle wavevector k remains unchanged (for- 
ward scattering), or it reverses its sign and changes by the 
amount of 2k (backscattering). As a result, backscatter- 
ing amplitude (in a single potential realization) vanishes 
in the Born approximation for 2k > q c , while in higher 
approximations, of order n, it vanishes for 2k > nq c . Cor- 
respondingly, in weak speckle disorder, LE is expected to 
exhibit a series of cascading drops at k — ^q c) so called 
"effective mobility edges". Thus, in order to assess the 
localization properties for 2k > q c , one should know LE 
beyond the Born approximation, which is a subject of 
this paper. 

LE in ID disorder was studied extensively for a variety 
of continuous and discrete models, which usually dealt 
with an uncorrelated disorder or a specific type of corre- 
lation (see Ref.[l3| and reference herein). For arbitrary 
correlation and in the weak disorder Born approximation, 
LE is known to be proportional to the disorder power 
spectrum E2>EH- In the recent years, correlations in 
ID disorder attracted much attention because they can 
induce unusual localization properties, such as existence 
of extended states 11311 . deviation from the single param- 
eter scaling (SPS) [14[ and appearance of the "effective 
mobility edge" when LE vanishes in Born approximation 
[HI, Eg. The latter means that exact localization length 
would exceed system size for sufficiently weak disorder, 
while higher orders are required to specify this regime 
quantitatively Ell- 
in this paper we report a systematic weak disorder ex- 
pansion for LE for two orders beyond Born approxima- 
tion, which is then applied to the laser speckle disorder 
E3- The analytical study is verified by numerical simu- 
lations, which also allow examination of regimes beyond 
weak disorder. The author of Ref-EB] performed a weak 
disorder expansion for a "generalized" LE. Let us stress 
that, although the generalized LE is easier to compute, it 
is the standard LE, studied here, that is of prime interest 



2 



in the localization problem. While these two quantities 
coincide in the lowest order [l6[, they differ in higher 
orders (see below). More generally, the equality holds 
under the assumption of SPS, and the differences found 
show how the latter is affected by disorder correlations. 
We consider continuous one-dimensional model 



+ fc 2 (l + 77(^ = 



(1) 



where rj(x) is the dimensionless disordered potential (for 
a quantum particle, rj(x) = —2mV(x)/h 2 k 2 is the ratio 
of the potential V(x) to the particle energy). We assume 
zero mean finite-range correlated disorder, specified by 
its (joint) cumulants (n — 2, 3, . . .) 



K"n (^1 , • ■ ■ , ^n — 1 ) 
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where R c is the correlation scale, the dimensionless func- 
tions r„ [xi, . . . , Xn-i) decay on the scale of unity and g 
is noise intensity. One chooses L °° r 2 (x)dx — 1, so that 
the limit R c — > yields Gaussian white noise with the 
two-point correlation (rj(x)ri(x')} = gS(x — x'). 
LE is defined as 



A= lim (2x)- 1 (ln(k' 2 \ip(x)\ 2 + \^'(x)\ 2 )) 



(3) 



where (..) denotes the disorder average, while the 
generalized LE considered in Ref.[l6| is given by 
lim ±- ln(fc 2 \ip {x)\ 2 + \ip' (x)\ 2 ). We calculate the LE 

x — too 

(JHJ) using the phase formalism relation 



(4) 



A = / zP st (z) dz, 



where z = ijj' /ip and P s t (z) is the stationary (i.e. the x — > 
oo limit) distribution of z. Introducing phase 9, defined 
by z — — k tan (0/2) and obeying "evolution" equation 



86 
dx 



2k + 2kcos 2 {9/2)rj(x), 



(5) 



LE is expressed in terms of the stationary distribution of 
the phase as 



X = - k J tan (5) P *t (0) d6. 



(6) 



Then, calculating the weak disorder expansion for P st (9) 
and substituting it into ©, yields the required expan- 
sion for the LE A. In the case in which i](x) in (O is 
a (5-correlated process, one can use standard technique 
[ll ] to write down a Fokker-Plank equation for P(9;x). 
When r/(x) is a correlated process, the method of "or- 
dered cumulants" (l9| can be used to obtain an approxi- 
mate master equation for P(0] x), given by a perturbative 
expansion in powers of /37 1 / 2 , where the dimensionless 
parameters 



7 = 2kR c and (3 = y/gk/8, 



(7) 



describe disorder correlation scale and strength respec- 
tively. Specific application of the method to the present 
problem involves many technical details, which will be 
given elsewhere [20| . Here we only outline the main steps 
of the derivation. Assuming small /3, but arbitrary 7, 
the following formal expansion of the master equation 
for P (6; x) is obtained 
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—P(e-x)= A0 + 2Y] (3 n K n P(9;x 

OX ^-^n = 2 



(8) 



K n are differential operators, whose definition involves 
integration over various ordered products of the operator 
Ax{x) = e xAo A l e~ xAo = -2kd e cos 2 (9 - kx), weighted 
with combinations of the cumulant functions T m < n (|2j). 
Here operators Aq — —2kdg and A-y — — 2kdg cos 2 (0/2) 
are related respectively to the deterministic and the 
stochastic terms on the right side of ([5]). 

Next, in the stationary limit, the master equation for 
P(9; x) reduces to an ordinary differential equation for 
Pst (9) — lim P(9; x) and one looks for a perturbative 

X — >oc 

solution 



P st (d)=Y,p n P n (6) 



(9) 



Let us note that this "double-stage" perturbation ap- 
proach is inconsistent in some cases [2l|,[22|, such as 
when P s t(9) is singular in the limit (3 — > 0. This, how- 
ever, does not occur in our case, since P s t{9) becomes 
uniform for (3 — > 0. Thus, substitution of the solution © 
into (J6j) yields the required LE expansion 



A = fc E >9 ^„( 7 ). 

' * n >2 



(10) 



The first four coefficients are given by 

A 2 = c , A 3 = c 6 , A 4 = c (ci - c 2 )+c 3 ci-(-2c4 + C5, (11) 
where Ci are the following functions of 7(= 2kR c ): 



i* 00 

Co = / dsT2 (s) cos 7s 
Jo 

c\ = J dsT2 (s) sin 7s = J 



T 2 (7) dc 
— , c 3 = - 7 —, (12) 



dq f 2 (q) 



, c 2 = 7- 

Ztt 7 — q 07 



ds\ I c?s 2 r 2 (si) r 2 (s 2 ) x 
Jo 

x [2 sin (7S1) - 7 (si - s 2 ) cos (7S1)] 



(13) 



1^ 



r 2 (7-9) , r 2 (7 - q) - r 2 (7) 

h T T~2 

q 2q 2 



p OO POO pOO 

c 5 = — 7 / dsi / ds 2 / ds 3 r 4 (si, s 2 , S3) x 
</ J 81 J 8% 

x (2 cos (7S3) + cos (7S1 - 7s 2 - 7S3)) , (14) 
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pOO POO 

CQ = -yJ 7 2ql ds 1 ds 2 T 3 (si,s 2 )sin7s 2 

JO JS! 

V27 f dq 



4tt 



r3 (g,7)- r 3 (q,i-q) 



(15) 



and r„ (gi, . . . , q n —i) is a Fourier transform of 
T n (xi, . . . , x„_i). The coefficients Cj for i = 5, 6 may be 
tagged " non-Gaussian" , since they depend on the higher 
(non-Gaussian) cumulants T3 and T4 only, and vanish 
for Gaussian disorder. Note that Ci = o,3,4.6 would van- 
ish for 7 above certain threshold, if T2 (<?) and T3 (91,(72) 
have finite support (same applies to C5, whose expression 
in terms of T± is not shown). Besides, if T2 (q) or its 
derivative are discontinuous at some point q c , then C4 (7) 
would diverge at 7 = q c and the perturbation theory 
would break down for this 7 (see below). 

As expected, A2, the lowest order coefficient in (fT0|) . co- 
incides with earlier results [l2j . Comparing our results for 
the standard LE to those of Ref. [l6| for the generalized 
LE (even orders only), one finds that the two coincide 
in the second and vary in the fourth order, though the 
difference exists only for non-Gaussian disorder. Namely, 
the corresponding non-Gaussian terms C5 (7) in A4 have 
different expressions, which at small 7 scale differently 
with 7 and have opposite signs. Thus, the SPS relations 
(see e.g. Ref. [HI) hold up to the order four in Gaussian, 
while are broken in non-Gaussian disorder. 

So far our results are quite general and pertain to an 
arbitrary random potential, the only condition being suf- 
ficiently fast decay of the disorder cumulants [H| . Now 
we specialize to the case of a laser speckle potentials, pro- 
duced by transmitting laser beam through a diffuser with 
a rectangular aperture Its intensity pair correla- 

tion function is To fx) = — sm 2 x , whose Fourier transform 
T 2 (q) vanishes for \q\ > 2 (F 2 (q) is related to the shape 
of the optical aperture, therefore for arbitrary but finite 
aperture it would always have finite support [7j). Assum- 
ing that speckle field (as opposed to intensity) is a com- 
plex Gaussian variable, which is true for sufficiently large 
diffuser Q , one concludes that any T n is expressed solely 
in terms of " irreducible" products of the two-point field 
correlators w (Xj — Xj). For example, the 3 rd cumulant is 

T3 (xi, x 2 ) = —2e (f ) 3 ^ 2 w (xi)w (xi — x 2 ) w (x 2 ), where 
w (x) — ^is^ is related to the Fourier transform of the 
"rectangle" function and e = ±1 is the sign of the disor- 
der distribution skewness, depending on cither "blue" or 
"red" laser detuning from the atomic transition. 

Substituting the explicit expressions for T n into the 
definitions of the coefficients Ci, one obtains the following 
LE expansion coefficients for the speckle disorder: 

A 2 = ^px(2- 7 ), A 4 = (A 4 G + A 4 WG )x(4-7), 



A, 



2s 



(7 - 2) In - 7 In | 



X(2~7), (16) 



where X G = Co (ci — C2) + c^Ci + 2c4 and Xf G = c 5 denote 
the "Gaussian" and the non-Gaussian parts of A4 and 



X (x) is the Heaviside step function. For 7 < 2 



A4 27 



4 7 - 3 7 " + ( 2- + 2 I In 



7 + 2 



+ [ ^--lOj ln^I+ 7 (7-4)ln| 



. (17) 



A-«f + l(2- 7) ln(2-7)(Hln^-l n 2 



while for 2 < 7 < 4 



\G _ 1 
A4 ~ 27 

\NG _ 7 
7T 



1 



47 - (7 - 47 



In 



7-2 



(18) 



2 L12 h m — 

76 2 



7(7-2) ln 2 7"_2 



where Li2 (z) is the dilogarithm and X^ G for 7 < 2 was 
calculated approximately, assuming (2 — 7) <C 1 and ne- 
glecting O (2 — 7) terms (exact calculation is too lengthy 
to carry out, while A^ G is important only close to or 
above 7 = 2). The coefficients Xi (7) are plotted in the 
insert of Fig. [2j As already noticed @, A 2 vanishes for 
7 = 2k R c > 2. Then, as can be expected from the general 
perturbation theory for scattering, A3 and A4 vanish for 
7 > 2 and 7 > 4 respectively. Thus, for 2 < 7 < 4, LE 
switches from quadratic to quartic leading dependence 
on (3 and, for weak disorder, undergoes a steep decrease 
at 7 = 2kR c = 2. Note that for speckle potential our 
perturbation theory breaks down at 7 = 2, where A G 
has a logarithmic divergence. As explained above, this is 
because T2 (q) = (2 — \q\) x (2 — \q\) has a discontinuous 
derivative at q = 2. 

Our analytical study was supplemented with numerical 
simulations of the discrete tight-binding model near the 
energy band edge, where it is a good approximation to 
continuous problem. LE was computed using the trans- 
mission matrix formalism [23| . LE, as a function of 7 at 
fixed 770 = \J (?7 2 ), the ratio between the disorder stan- 
dard deviation and the particle energy, is shown in Fig. 
Q]for r]o = 0.08 and r] = 0.15 (the disorder strength pa- 
rameter (3, rewritten as j3 = j^/T^yrjQ, changes along with 
7 in this parametrization). LE, for each r]o computed for 
both signs of e, is significantly larger for e = +1, which 
demonstrates strong effect of the non-Gaussian character 
of the speckle disorder. While agreement between the an- 
alytical and the numerical results is acceptable for 7 < 2, 
it appears to be very poor for 7 > 2, which is because 
the disorder is not weak enough. To clarify this point, 
numerical LE (A^™"'), computed for different values of 
(3 and e at fixed 7 = 2.2 and divided by the analytical 
A = fc/3 4 A4, was fitted to a 4 th degree polynomial of e/3, 
as appears in Fig. [5] together with the fit equation (re- 
call that for 7 > 2 the perturbative expansion (fT0|) for 
A starts with the fourth order term). The free constant 
of the fit is \( 4 Num '> y_\ 4) anc [ its value 1.003 indicates ex- 
cellent agreement between the perturbation theory and 
numerics. On the other hand, the fit equation in Fig. 
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FIG. 1: (color online) Lyapunov exponent A for r/o — 0.08 
(main panel) and 770 = 0.15 (insert) as a function of the di- 
mensionless disorder correlation parameter 7 = 2kR c . Cir- 
cles and squares show numerical data for positive (e = +1) 
and negative (s = — 1) skewness of disorder, while solid and 
dashed lines give corresponding 4 th order expansion. Curve's 
spike at 7 = 2 reflects the logarithmic divergence of Aj. 
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FIG. 2: (color online) LE calculated at fixed 7 = 2.2 as a 
function of disorder strength and "sign" e/3. Circles represent 
numerical LE devided by fc/3 4 A4. Line and equation show 
the 4 th degree fit over all the points except the left- and the 
right-most ones. Insert: expansion coefficients Ai (7). 



[5] shows that the 4 th order approximation is acceptable 
only for very weak disorder (/3 <C 0.1). Repeating similar 
test for 7 a 2, we conclude that the logarithmic peak of 
A4 (Fig. [21 insert) is reproduced in the numerics up to 
the smearing effect of the finite system size L (this can 
be accounted for by smoothing Ai (7) over A7 ~ -^). 

Note the very rapid growth of the higher order coeffi- 
cients of the fit in Fig. [21 corresponding to higher orders 
in LE expansion. As suggested by the structure of the 
perturbation theory, it is related to the non-Gaussian 
character of the exponential distribution of the speckle 
intensity, whose cumulants grow factorially. This ex- 
plains both the failure of the 4 th order approximation 
and the strong dependence of LE on the sign of the dis- 
order skewness e found at moderately weak disorder (Fig. 
[T]). The perturbative expansion up to order /3 4 accounts 
for the effect of e only for 7 < 2, since its only odd term 
A3 vanishes for 7 > 2. Therefore, analytical curves for 
e = +1 and e = — 1 coincide for 7 > 2 (Fig. [T]). 

Finally, we address persistence of the effective mobil- 
ity edge at 7 = 2 as a function of the disorder strength. 
This question is important for the recent experiments on 
the BEC expansion @ (with e = +1), because, if dis- 
tribution of the atomic fc's in the "exploded" condensate 
stretches beyond the speckle frequency cutoff 2/R c , then 
the deepness of the LE drop at 7 = 2 becomes crucial 
for predicting algebraic versus exponential decay of den- 
sity profiles [3]. Our computations in Fig. [T]show that 
this drop depends strongly on the disorder strength and 
becomes of effectively one order of magnitude or less for 
rjo > 0.1 (relevant to the experimental regimes Q). In 
addition, note the difference in LE for "blue" (e = +1) 
and "red" (e = —1) detuning: for negative detuning LE 
is smaller, but its relative variation across 7 = 2 is larger. 



(Though our parametrization is not natural for this ex- 
periment - we fixed k with disorder amplitude and varied 
i? c , instead of varying k with the rest fixed, - the conclu- 
sions remain valid, since mainly the behavior near 7 = 2 
is concerned.) 

In conclusion, we derived general expansion for the 
Lyapunov exponent (LE) in ID correlated disorder two 
orders beyond the Born approximation. Comparing it 
with that for the generalized LE [lj| shows that single 
parameter scaling is broken perturbatively in the fourth 
order for non- Gaussian disorders (with finite moments). 
Applying this expansion to speckle disorder with Fourier 
spectrum cutoff q c = we find that the leading order 
dependence of LE on the disorder strength crosses from 
quadratic for kR c < 1 to quartic for 1 < kR c < 2, as ex- 
pected from the standard QM perturbation theory. For 
very weak disorder, this results in large and steep drop 
of LE across kR c — 1. For larger but still weak disorder, 
this drop moderates, while LE becomes very sensitive to 
the skewness of the disorder distribution (for kR c > 1). 
This is because of strongly non-Gaussian distribution of 
speckle intensity, whose cumulants grow factorially fast. 
Physically, this reflects non-perturbative contribution of 
rare but large potential peaks (dips) of the typical width 
R c to the scattering process. 
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